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Phase Matrix Induced Symmetries for Multiple Scattering Using 
the Matrix Operator Method 

By Stephen J. Hitzf elder and George W. Kattawar 
Department of Physics, Texas A&M University, College Station, Texas 77843 

I 

ABSTRACT 

Entirely rigorous proofs of the symmetries induced by the phase matrix into 
the reflection and transmission operators used in the matrix operator theory are 
given. Results are obtained for multiple scattering in both homogeneous and 
inhomogeneous atmospheres. These results will be extremely useful to researchers 
using the method in that large savings in computer time and storage can be 


achieved. 
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Introduction 

Symmetry relations play an important part in methods of solving radiative 
transfer problems. They can be used to reduce the amount of calculation and core 

storage and can also be used for giving physical insight into the results 
obtained. Principles of invariance played an important part in Chandrasekharas 
(1960) solution to the Rayleigh scattering problem. Other work involving S3nrjnetry 
relations include Sekera (1966) and Busbridge (1960). Recently, Hovenier (1969) 
proved that certain symmetry properties of the phase matrix induced equivalent 
symmetries in the scattering and transmission operators. The method used by 
Hovenier for proving these was cumbersome and only presented in the case of a 
homogeneous atmosphere. Kattawar (1973) derived the symmetry relations induced in 
the reflection and transmission operators in the matrix operator (MO) method from 
symmetries of a phase function.^ Tanaka (1971) proved a set of symmetry relations 
for the diffuse scattering and transmission operators. These results were in 
agreement with Hovenier *s. It is the purpose of this paper to extend MO to include 
polarization and to derive symmetry relations induced in the reflection and 
transmission operators using the MO method. The proofs will be rigorous for both 
the homogeneous and inhomogeneous atmqspheres; 

Matrix Operator Theory with Polarization 
MO has been developed over the years by contributions from many sources. A 
brief history of the development of the method and most of the important equations 
are given by Plass, Kattawar, and Catchings (1973). This work will be referred to 
as PKC. We present the equations in the continuum. The I, Q, U, V Stokes vector 
representation (Chandrasekhar, 1960) will be used in this paper. Therefore, we 


define: 
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I (T;y,(f>) = 


Q(x;y,(t)) 0<y<l (la) 

U(T;y,(J.) 

_V(T;y,(|)) ^ 

I(T,-y,<|i) 

Q(x;-y,(|)) 0<y<l (lb) 

U(x;-y,(J)) 

V(x;-y,(f>) 


as being the radiance in the direction of increasing and decreasing optical depth 
respectively. In Eqs. (1), y = cos 6, 0 being the polar angle measured from the 
zenith, and being the azimuthal angle where we use the external source as the 
reference for azimuth. We consider a plane parallel atmosphere with boundaries 
X and y (x<y) . Following PKC we have: 

^ ^ ^ ^ 1 ^ I ^ It 

I (y;y,<?>) = t(x,y;y ,<>;y ,<>) * I'^(x;y,(;)) + r (y ,x;y , ;y ,<J>) * I“(y;y,<{>) (2a) 

I (x;y,<>) = r(x,y;y,(J);y ,1) r l'*’(x;y,4)) + t(y ,x;y ,<J.;y ,«5) * I (y;y,(J) (2b) 

r 

In Eqs. (2) and in the sequel to follow, the • implies the following operation; 

t(x,y;y,(|);y ,{f>) • l'^(x;y,(J)) ^ /^dy/^^d^ t (x,y ;y ;y ,|) l'*'(x;y ,<!)) (3) 

^ . ■“ o o - 

+ - 

Subtracting I (x;y,<{)) and I (y;y,(J)) from Eqs, (2a) and (2b) respectively and passing 
to the limit x-^, the following differential equations are obtained: 

= -i*^(y;y»<l’;y »4>) * j"*'(y;y,<i>) + r'^(y;yi<l>;y>4>) ’ i (y;y,4>) (^a) 

” 3y = r (y;iJ'»4>;y>'i’) * I (y;y»«l>)*-- r (y;y ,4>;y,(f>) • I~(y;y,(f)) (4b) 
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where the gamma operators are given by; 


I = Lim(5(y-y)(5((f>-(f))E - t(x,y;p,4>;iJ,<{)))/(y-x) , 

X")-y 


(5a) 


r (y;p,<f);u,(|)) = Lim r(y,x;y,(|>;u,<|)) /Cy-x) , 
x->y 


(5b) 


I (y;y»<J>;y,<i>) = Lim r(x,V;y,<{);y,i) /(y-x) » 
x-^y 


(5c) 


r (y;y,<{);y ,^) = Lim ( 6 (y~y) 6 (<>-(}))£ - t (y ,x;y ,(^>;y ,(?>) ) /(y-x) . 
x*>y 


(5d) 


In Eqs. (5), E is a 4x4 unit matrix and 6(y-y) is the Dirac delta function. The 
gaimna operators, defined in Eqs. (5) are the generators of the r and t operators. 
Eqs. (4) form the link between the equation of transfer and the interaction 
principle of MO (Eqs. (2)). 

From Chandrasekhar (1960) the equation of transfer in a plane parallel 
scattering atmosphere, taking polarization into account, is 


d T 


+1 


2tt 


y + I(T,y,<{)) = w(t) / / p(Tjy,(});y ,(})) I(x ;y ,<{i)dyd(f) (6) 

-1 o " 


- f f 

where m(x) is the single scattering albedo, and p(x;y ,4> ;y ,4>) is the phase matrix 

where the element p^^ is normalized to unity. Using Eqs. (1) for 
+ — 

I ('c>y»<l>) snd I (x;y,<j)) we have from Eq. (6) after some simple algebra: 

d I^(x*u (b) ^ ’ * I I + f f • * 

— - ’ * = / / “(1/y) (6(y-y)6(({)-4))E - m(x)p(x;p,4);p,j>)) I (x;y ,i)dyd(j) 


dx 

o o 

+ /^^(co(x)/u)p(x;y,c});“p,4)) l“(x;y,4>) dy* d(|)^ 

o o 

_ I (t.M.lf) ^ dp ' dlfl ' 


(7a) 


dx 


o o 


.1 ,2tt 


t f t t 


+ / / -(1/y) (6(y-y)6(.i>-(|i)E ~ o)(x)g(x;-y ,({);-y ,4>)) l'(x;y ,<J>)dyd<J) 

o o 


(7b) 
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Comparing Eqns. (7a) with (4a) and (7b) with (4b) it is easy to make the following 
identification of terms: 


r (y;v,i^;y,^) - (6(p-ii)6(^-(^))E ~ m(y)p(y ,<j>) ) /p , (8a) 

£ (y == («S(p-y)6((f)-<}))E - 0 >(y)p(y;-y,(f);-p,<j))) /y , (8b) 

£ (y;p»^;y ,<l>) = t^(y) .p(y;v,<|!;-u,^)/y , (8c) 

^ ® t I 

£ (y;y,<l»;y,<t>) = w(y) p(y;-y,4>;y ,<j:)/y . (8d) 


Thus, the generators for the r and t operators are related to the single scattering 
albedo and the phase matrix. 


PKC give the equations that combine' the r and t operators for two arbitrary 
layers to yield the r and t operators for the combined layers. Using these equations 
to add an infinitesimal layer to the top and to the bottom and then passing to the 


limit as the small layer goes to zero, the differential equations for the r and t . 
operators are obtained (Kattawar, 1973). The equations are as follows: 


3r(y,x;y,(J.;y ,<j) ) 

— o o _T~ 


,+4* 


9y 


- £ ” £ (y;y,<});y ,<>) ’ r(y,x;y , 4 );y^, 4 )^) 


^ I I 

- r(y,x;y,(|>;y,(j)) • T (yiP ,<f-;y^,<f)^) 


” f- nut* It 

+ r(y,x;y,<f>;y,({)) -T (y;y ,<>;y ,<J>) • r (y,x;y ,<f> ;y^, (j)^) ^ 


3r(y,x;p,(j>;y^,<|,^) 


H — II II I t 

• _ * / 


9x 


t (x ,y ; y , (ji ; y , (^) p (x;y ,<j>;y ,(j>) b(y ,x;y y^, (|) q) ^ 


(9a) 

(9b) 


3r(x,y;v,<j>;y^,({.^) 

3y 


nut I 


I t 


= t(y,x;y,<J);y,0 * T (y ;y ,(|)) * t(x,y;y ,<}>;y^,<|)^) , 


(9c) 
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^ = r (x;y,(});v^,(i>^) - r (x;y ,(j) ;y ,<{>) ’ r(x,y;y,({);y ,<J> ) 


' ' "Ht I I 

- r(x,y;y,<|);y,0 * T (x;y ,({);y^,(})^) 


tt II t 1 


+ r(x,y;y ,<^;y,(j>) * V (x;y ,(J) ;y ,4>) * r(x,y;y ,<j>;y^,(f>^) , 


(9d) 


3t(y,x;y,iJ>;y^,(J.^) 

9y 


* ^ f T I I tl IT 

■t(y,x;y,(|);u,(})) • (p (y;y ,4>;y^,'i)^) - r (y;y,4>;y,4>) 


It ft 

* r(y,x;y,(});y^,(j)^) ) , 


(9e) 


8t(y,x;y,(J>;y^,(J)^) 

9x 


It II II _j II II I f 

(I' (y»y»^>U 30 ) “ ^ j y » y j ^ > y j r (xjy,(ji jy 


t I 

t(y,x;y,<^;y^,(j)^) , 


(9f) 


3t(x,y;y,(|.;y^,(})^) 

3y 


•([ (y;y ,<J>;y ,<>) - r(y,x;y’,<|);y ,(J)) * T (y ;y ,i|);y ,<j)) ) 


’ t(x,y;y ,4>;y^,(J.^) , 


(9g) 


8t(x,y;y,(^);y^.,(i)^) 

3x 


I t -H* If -f— I I ” ♦* 

■t(x,y;y,({);y,(J>) • (F (x;y ,<> ;y^,(|)^) - F (x;y ,tf) ;y ,4>) 


r(x,y;y ,(f;y^,(J)^)) . 


(9h) 


The boundary conditions for r and t are r(x,x) = r(y,y) = 0, and t(x,x) = t(y,y) = E. 
,Eqs. (9) may easily be sho\m to be equivalent to the differential equations derived 
by Chandrasekhar for the diffuse scattering and transmission matrices. The obvious 
difference of course is that the t operator of MO accounts for the direct and diffuse 
radiation whereas Chandrasekhar's transmission operator accounts only for the diffuse 
radiation. Thus, the differential equations for the t operators ((9e) - (9h)) are 
simpler than Chandrasekhar's corresponding equations. Also, the equations that 
couple r and t ((9b) and (9c)) are also simpler than the corresponding equations for 
the diffuse operators. 
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Symmetry Properties of the Phase Matrix and Gamma Operators 


The syrametry properties of the phase matrix have been used extensively in the 
past. Van de Hulst (1957) used symmetry alone to reduce the number of independent 
elements in the phase matrix. Sekera (1966) derived a set of symmetry relations 
for the phase matrix. Hovenier (1969) developed the symmetry relations for the 
phase matrix both from a space- t inie ’ symmetry basis and from a mathematical basis. 
Kattawar, Hitzfelder, and Binstock (1973) explicitly demonstrated these 
symmetries for the phase matrix of a polydispersion of spherical particles- Their 
relations were derived directly from Mie theory. These symmetry properties are 
reproduced in the notation of this paper for clarity. 


f 1 

p(p,<f>;d,'l>) = Q 

t t 

t I 

p(-U,(J>;-U,(J)) = 

I t 

p(y,<>;y ,(|>) = P 

I . t 

. p(y,<f>;y,<f>) = p 

( 

I X t 

I t 

p(-y,4>;-u,<l») = 

where 


rp f f 

p (u,(|);y,4>) Q , 

(lOa) 

ITI t f 

P p (y,4;;y,<{») p > 

(10b) 

1 1 

p(y,4>;y >«t>) > 

(10c) 

1 1 

Q p(y ,<r>;y ,<J») Q p » 

(lOd) 

r 

•Ti f I 

p (y,<f^;y,4>) p » 

(lOe) 

r|i t t 

Q p (y,<i>;y ,<})) Q , 

(lOf) 

P Q p(y,i|);y ,(J)) Q P , 

(lOg) 


“ 1 

0 

0 

0“ 


" 1 

0 

o 

o 

» 

0 

1 

0 

0 


0 

1 

0 0 

0 

0 

-1 

0 


0 

0 

1 -0 

_ 0 

0 

0 

1 


_^0 

0 

0 -1 
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and the superscript T implies a matrix transpose. Eqs* (10) form the complete set 

of possible transformations that can be performed on p (y i ;p ,<J>) . As Hovenier points 

out in his derivation only three of the seven relations are independent. The other 

four may be derived in terms of a basis set. This idea can be deduced very simply 

if one realizes that the seven transformations on the phase matrix coupled with 

the identity transformation, which leaves the phase matrix unaltered and denoted 

by E, forms an abelian group of order eight. The group multiplication table is 

given in Table I, A close inspection of Table 1 shows that there are ten possible 

subgroups of order four. These subgroups can be found as follows. Let a and 3 be 

any two distinct elements of Table 1. Then a, 3, end a3 will form a proper subgroup 

of order four when the identity element is included. The reason for this is clear 

2 

since every element of the group is of second order, i.e., a = E where a is any 
element of the group. Therefore, a(a3) = 3 and 3(«3) = ct. This will leave tw^enty 
five triplets of elements wrhich can be used as generators for the entire group. 
Hovenier only considered the triplets (d,e,g) and (b,c,d). From the space-time 
symmetry point of view, (b,c,d) form the basic set of relations. 

The symmetry properties of tlie gamma operators are now derived from Eqs. (10) 
and Eqs. (8), Transposing Eq . (8a) we have: 


= (5(p-|!i)S«.4)E - u(y)p’’^(y ;y ,0 /n 


Pre.and post multiplying by Q and 


' * * 

notine that fifu-u^ = 5fu-u) 


r’’^(y;y ,4>) = Q (r^(y;y Q 


Similarly from Eq, (10a) and Eqs. (8b), (8c), and (8d) 

r““(y;y,<{>;y,(J)) = Q (r (y;y ,4>))*^ Q (lib) 

= Q (f”^(y;v,<l>;u,<l>))'^ Q (iic) 


Similar operations can be performed on Eqs. (8) taking into account Eqs. (10) and 
the following relations can be derived: 
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I I 

r = 

j<l>) = 

-+ ' ' 

r (y;y,<f>;y,<[)) = 

r"*~^(y;y,<{>;y,'+>) = 

+- * ’ * 

r (y;u><{>;y ,<i>) = 

t I 

' . !”■ (y;y,(^;y,4>) = 

+- ’ * 

r (y;y,(+.;y,(J)) = 

r (y;y,<i>;y,({») = 

++ * * 

r (y;y,<ij;y,(J') = 

I 1 

(y 7 y » ^ > y > *t’) ~ 

+- ' * 

r (y;y,<j»;y,(}>) = 

4 

I t 

r (y»y>^N»y»^) ~ 

+- ' * 

r Cy;y,(J);y,(j>) = 

-I- * ' 

r Cy;y,cJ);y,(})) = 

++ ' ’ 

r (y;y,(}>;y ,(j)) = 


I (r‘^(y;y,'j':y,«)''^ p . 

P (r‘^(y;y,'l':w,«)^ p, . 

P (P~'’(y;u,ii>;y,'{i))^ p > 

f I 

r (y;y,<l’;y,4>) » 

-+ * * 

r (y;y,'J>;y ,<>) » 

P Q r‘^(y;y,'t>;w,« Q p • 

P Q r (yiy ,4>;y ,4>) Q P , 

+- ' ’ 

PQr (y;y,tJ);y,<{))QP, 
P Q r"*'(y;y.«;v.«) Q p . 

P(r (y ;y ,(|);y ,ip)) p » 
p (r“"(y;y,<);i,i))'^ P , 
P(r~^(y;y,(t>;y,<))^'^ P , 
Q(r~^(y;y ,<>;y,(J>))'^ Q , 
Q(r"*^(y;y,<i>;y Q , 
Q(r^(y;y,<f>;y Q , 

P Q r (y;y,<i>;y ,4.) Q P , 

P Q ^”^(y;y,<^;y,<^>) Q p . 
~ ~ ^ ^ 


(12a) 

(12b) 

(12c) 

(13a) 

(13b) 

(14a) 

(14b) 

(14c) 

(14d) 

(15a) 

(15b) 

(15c) 

(16a) 

(16b) 

(16c) 

(17a) 

(17b) 


The reader may easily verify that Eqs. (12) to (17) were derived from Eqs. (10b) to 
(lOg) respectively. The interdependence of Eqs. (10) apply to Eqs. (11) to (17) as 
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well. 


Symmetry Properties of the r and t Operators 
The symmetry properties of the gamma operators are a necessary step in 
deriving the properties of the r and t operators. These properties will now be 
derived in a straight forward manner for both the homogenous case and the 
inhomogeneous case. Transposing Eqj (9a) and pre and post multiplying by Q 
one obtains: 

Q Q(r"^(y ;y Q 

” Q r*^(y,x;y,(});y^,<ji^)Q * Q(r'^(y;y,4>;p Q 
" Q(r Q * Q r'^(y,x;y,^;y ,<J») Q 

+ Q Q • Q(r ' (y Q * Q r^(y ,x;y ;y ,4>) Q (IS; 

~1 

where Q = Q has been used. Using Eqs. (11) and letting u(y ,x;y ;y^,4>^) = 

T 

Q r (y ,x;y,(J);p^,(j)^) Q one obtains: 

g f. * • -j. { * ♦ 

u(y,x;y,4>;y ,<}) ) = T (y p ,4>) - u(y,x;y,(|);y ,<}) ) * r (y ,y ,(J>;y ,({)) 

Oy ^ OO^ UO .w OU ,-^ 

t t It 

- r (y;y^,<j>^;y,4>) * u(y,x;y ,<});y ,(J>) 

I t _| I I II IT 11 II 

+ u(y,x;y,(});y^,(})^) ‘ r "(y ;y ,(f);y ,<}>) * u(y ,x;y ,«J>;y ,<|)) (19) 

Comparing the above equation with Eq. (9d) it is easy to. see that if the gamma 
operators are independent of optical depth, that is the atmosphere is homogeneous, 
then u(y ,x;y ,4> ;y^,<i>^) satisfies the same differential equation as r(x,y ;y^ ,4>^,y (f;.) 
and has the Scune boundary conditions. Thus, the solution is unique and we get 




( 20 ) 
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Similar operations can be performed on Eqs. (9). Another example is presented to 
demonstrate the derivation of a symmetry relation in an Inhomogeneous atmosphere. 
Transposing Eq. (9a) and pre and post multiplying by P one obtains: 

^ P = P (r*^(y;y P 

- P (r (y;y' :P q.<!>q))^p ■ p r^(y,x;tj,ii>;tj' , 4 .') p 

+ P r’'^(y,x;p',4.';u^,<j,^) p • p (r‘^(y,v",((.":y’ , 4 .'))’^ P 

• P r'^(y,x;u,4.;y",4>") P (21) 

where P = P ^ has been used. Using Eqs. (12) and letting w(y ,x;y ,<}>;y^,(|)^) = 

T 

P r (y ,x;y., (J>;y^,(|>^) P one obtains: 

■“ w(y,x;y,<>;y^,^^) - (y;y^,<{>^;y ,<>) - w(y ,x:y * ,(J> ' ;y^ ,^^) * T (y J y * ,<J> ' ;y ,4>) 

- r (y;y^,(|>^;y' ,(|>*) ’ w(y ,x;y ,^>;y ’ ,^*) 

+ w(y,x;y ',({)' ;y^,(^)^) • r*^(y;y \(j) ’ ;y",(J)") * w(y,x;y,<|);y",(J)'‘) (22) 


Comparing the above equation with Eq , (9a) it is easy to see that 


r(y,x;y^,<^^;y,(Ji) 


P r (y,x;y ,0;y^,(|)^) P 


(23) 


The above relation will hold in any plane parallel atmosphere since the gamma operators 
are evaluated at the same level in both differential equations. It is a very simple 
task. to perform similar operations on Eqs. (9) and using Eqs. (11) to (17) to derive 
the complete set of symmetry relations for the r and t operators. As in the two 

examples given above, it is easy to determine whether the relation is valid in an 

/ 

■inhomogeneous atmosphere or only in a homogeneous atmosphere. The complete set of 
relations is given below. In this list, an i(h) means that the relation holds in an 
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inhomo gene ous (homogeneous) atmosphere. 



T 

= Q r (y,x;y,4);y^,(f)^) 


h 

(24a) 


= Q t"^(y,x;y,4);y^,(}.^) 

Q . 

h 

(24b) 


T 

= Q t (x,y;y,4>;y^,4>^) 

Q . 

h 

(24c) 

« 

r(y,x;y^,4>^;v,({)) 

T 

= P r (y,x;y,4>;y^,(j)^) 

P , 

i 

(24d) 

r(x,y;p^,<j>^;y,<|)) 

T 

= P r (x,y;y,4.;y^,(j)^) 

1 , 

i 

(24e) 

l(x,y;p^,(f>^;y ,4>) 

T 

= P t Cy,x;y,(f);y^,4>^) 

P , 

i 

(24f) 

r(y,x;vi , 

= r(x,y;y,4>;y^,<>^) , 


h 

(24g) 

t(y,x;y,({)^;y^,<J)) 

= t(x,y;y,4»;y^,4)^) , 


h 

(24h) 

r(y,x;y,({)^;y^,(J>) 

= P Q r(y,x;y,<>;y^,<|)^) 

Q P , 

i 

(24i) 

r(x,y;y,(l)^;y^,c;)) 

= P Q r(x,y;y,({);y^,4)^) 

Q P , 

i 

(24j) 

t(x,y;y,i^)^;y^.,4)) 

« 

= P Q t(x,y;y,4>;y^»<J’Q) 

Q P , 

i 

(24k) 

t(y,x;y,(})^;y^,<|)) 

i 

= P Q t(y ,x;y,4>;y^,4.^) 

Q P , 

1 

(24I-) 

r(x,y;y^,ij);y,^^) 

T ' 

= P r Cy,x;y,4);y^,4)^) 

P , 

h 

(24m) 

t(y,x;y^,(J.;y ,4>^) 

T 

= P t (y,x;y,4>;y^,4)^) 

P , 

h 

(24n) 

h (x ,y >y^j4*>yj 

T 

= P t (x »y Jy , 4> 5 y^ , 4»^) 

P , 

h 

(24o) 

r(y,x;y^,({>;y ,4)^) 

T 

= Q r (y,x;y,4);y^,4.^) 

Q * 

i 

(24p) 

r(x,y;y^,4;;y,4>^) 

T 

= Q r (x,y;y,(|);y^,4.^) 

Q > 

i 

(24q) 
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= Q t^(y ,x;y,4>;y^,(|i^) Q , 1 (24r) 

r(x,y;y,<|);y^,({>^) = P Q r (y ,x;y ;y^,4>^) Q P , h (24s) 

t(x,y;ii,4>;y^,(J)^) = P Q t(y ,x;y ,(f);y^,(i.^) Q P . h (24t) 

The above set of symmetry relations has been rigorously proved for both the 
homogeneous and the inhomogeneous cases. They are identical to the relations derived 
by Hovenier (1969) and can be shown to be identical to those derived by Tanaka 
(1971) for the diffuse scattering and transmission matrices in another Stokes vector 
representation. The method of proof is straightf or^-jard and yields physical insight 
into the symmetry relations. 


M 



Conclusions 


It has been sho^^m that the differential equations for the reflection and 
transmission operators used in matrix operator theory can be used to derive 
sjrmmetry relations induced in them by the phase matrix. This technique is entirely 
rigorous and straightfor^v^ard. One of the most powerful results of it is that the 
results for multiple scattering in inhomogeneous media are as easy to establish as 
those for homogeneous media. These symmetries, v;hen properly used, can result in 
vast savings in computer time and storage. This is essential if one desires to 
use the method for multiple scattering in a realistic atmosphere where inhomo- 
geneities are encountered. 
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Table 1. Group multiplication table for the complete set of 
transformations on the phase matrix. 
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